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METHOD AND APPARATUS FOR MODELING 
COSMIC RAY EFFECTS ON MICROELECTRONICS 

BACKGROUND OF THE INVENTION 

Field of the Invention 

The present invention relates generally to a method and apparatus for modeling the effects 
of cosmic rays on microelectronics in earth orbit, and more particularly to an improved method and 
apparatus for modeling these effects, these improvements reflecting a simpler, easier way for the user 
to operate the model over the internet, and more accurate modeling of these effects. 

Description of the Related Art 

For electronic components onboard satellites in earth orbit, exposure to cosmic rays 
represents a serious risk, due to the capacity of cosmic rays to induce single event effects (SEE) in 
these components. See generally Sherra E. Kerns, Transient-Ionization and Single-Event Phenom- 
ena, in Ionizing Radiation Effects in MOS Devices and Circuits 485-9 1 (John Wiley & Sons, 
Inc., T.P. Ma et al. eds., 1989). In brief, SEEs occur when an energetic particle changes a particular 
device in an integrated circuit, thereby causing an error. To date, the only effective methods for 
hardening these circuits against SEEs have been shielding, redundancy, and error detection and 
correction (ED AC). 
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Unfortunately, these measures are not always effective, and the failure of a single component 
may lead to the total loss of a multi-million dollar satellite. With the potential loss associated with 
this risk so high it is desired to have a way to accurately predict the magnitude of this risk, as well 
as a way of predicting how successful ameliorative efforts are likely to be. 

CREME (short for Cosmic Ray Effects on Microelectronics) was a software package 
developed by the Naval Research Laboratory in 1981 for modeling how a given electronic chip on 
a given satellite with a given orbit and amount of shielding would hold up against cosmic ray 
bombardment. See Cosmic Ray Effects on Microelectronics, Part I: The Near-Earth Particle 
Environment, Adams et al., NRL Memorandum Report 4506; Cosmic Ray Effects on Microelectron- 
ics Part II: The Geomagnetic Cutoff Effects, Adams et al., NRL Memorandum Report 5099; Cosmic 
Ray Effects on Microelectronics, Part IV, Adams, NRL Memorandum Report 5901, each incorpo- 
rated herein by reference, in their entireties, for all purposes. CREME had several shortcomings. 
It has been discovered that many of the predictions of CREME were inaccurate. Particular 
shortcomings of the original CREME software included its inaccurate modeling of the transmission 
of cosmic rays through earth's magnetosphere, and inaccurate modeling of the flux of heavy ions 
associated with solar flares. 

Moreover, the implementation of this program was less than optimal, having a difficult user 
interface, and requiring each user to install, maintain, and run the software. 
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NOTE ON NOMENCLATURE USED HEREIN 
Some portions of the following disclosure describe method elements such as determining, 
selecting, dividing, processing, computing, calculating, numerically integrating, applying a function, 
displaying, and the like, and structure elements such as tables, memories, internet connections, and 
the like. These descriptions are the means used by persons of ordinary skill in the art of data 
processing to most effectively convey the substance of their work to other skilled artisans. Such 
methods and structures are intended to describe methods and structures for carrying out a set of steps 
on at least one programmed digital computer to reach a desired result. Thus, each of these steps 
requires a physical manipulation of concrete quantities, generally in the form of electrical, optical, 
and magnetic signals capable of being stored, retrieved, combined, and otherwise manipulated in 
such a programmed digital computer. 

Accordingly, unless indicated otherwise, skilled artisans will recognize that as used herein, 
terms such as determining, selecting, dividing, processing, computing, calculating, numerically 
integrating, applying a function, displaying, and the like, refer to the operations of a programmed 
digital computer system, or similar electronic computing device, that manipulates and transforms 
data represented as physical quantities within the computer. 

SUMMARY OF THE INVENTION 
Accordingly, it is an object of this invention to improve the modeling of the effects of cosmic 
rays on microelectronics. 

It is a further object of this invention to improve the modeling of solar heavy ion flux. 
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It is a further object of this invention to improve the modeling of geomagnetic transmission. 

It is a further object of this invention to improve the user interface, by implementing an 
internet-based interface, allowing users to access a world wide web site connected to a server where 
the software is installed and maintained. 

These and additional objects of the invention are accomplished by the structures and 
processes hereinafter described. 

An aspect of the present invention is a method and apparatus for computing a geomagnetic 
transmission function. This apparatus includes a programmed digital computer running modeling 
software for modeling the transmission of cosmic ray particles through the magnetosphere. The 
software includes a model representing a solution to the Lorentz equation in a magnetic field given 
by: 

•B = B IGRF (r,t , ) + B TSYG (Kp,r,f) 
where B is the earth's magnetic field, where B IGRF (r,f) is the International Geomagnetic Reference 
Field model of earth's magnetic field, a standard internationally recognized representation of earth's 
magnetic field, and where B TSYG (Kp,r,t*) is the model of earth's magnetic field published by 
Tsyganenko in 1989, as modified by the inventors. B IGRF (r,f) and B TSYG (Kp,r,t') are discussed in 
further detail infra. 

Another aspect of the present invention is a method and apparatus for computing a flux of 
particles at the outer surface of a satellite comprising, inter alia, an improved method and apparatus 
for computing a flux of solar heavy ions. This apparatus includes a programmed digital computer 
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running modeling software for modeling the flux of cosmic ray particles through the outer surface 
of a satellite. 

Another aspect of the invention is a method and apparatus comprising a programmed digital 
computer running modeling software for modeling the effect of cosmic rays on microelectronics, 
where this software embodies at least one of the two foregoing aspects of the invention. 

Another aspect of the invention is a preferred embodiment of a method and apparatus 
comprising a programmed digital computer running modeling software for modeling the effect of 
cosmic rays on microelectronics, where this software embodies at least one of the two foregoing 
aspects of the invention, where this preferred embodiment is connected to a network, typically the 
internet, to permit remote users to use the invention. 

, BRIEF DESCRIPTION OF THE DRAWINGS 

A more complete appreciation of the invention will be obtained readily by reference to the 
following Description of the Preferred Embodiments and the accompanying drawings in which like 
numerals in different figures represent the same structures or elements, wherein: 

FIG. 1 is a partial flowchart for a preferred embodiment of the invention. 

FIG. 2 is an elevation of the earth, showing the path of a cosmic ray particle in the 
magnetosphere. 

FIG. 3 is a flowchart for a method for calculating a geomagnetic transmission function for 
cosmic rays according to the invention. 
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FIG. 4 is a plot of a geomagnetic transmission function for a low inclination space shuttle 
orbit under quiet magnetospheric conditions. 

FIG. 5 is a plot of a geomagnetic transmission function for a lo inclination space shuttle orbit 
under stormy magnetospheric conditions. 

FIG. 6 is a plot of a geomagnetic transmission function for a space station orbit under quiet 
magnetospheric conditions. 

FIG. 7 is a plot of a geomagnetic transmission function for a space station orbit under stormy 
magnetospheric conditions. 

FIGS. 8 and 9 together are a flowchart for a method for modeling solar heavy ion flux 
according to the present invention. 

FIG. 10 shows a main menu for a user to operate the preferred CREME96 software through 
the internet. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
The following are each incorporated herein by reference, in their entireties, for all purposes: 

(A) "CREME96: A Revision of the Cosmic Ray Effects on Micro-Electronics Code", 
Transactions on Nuclear Science 44(6) 2150-60, Tylka et al. (1997); 

(B) "A Magnetosphereic Magnetic Field Model with a Warped Tail Current Sheet", Planet. 
Space Sci. 37(1) 5-20, Tsyganenko (1989). 
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Overview 

In a preferred embodiment, the invention comprises several computer software modules 
running on a powerful server computer. Although this computer may be operated in a standalone 
mode, it is preferred to have this server connected to a network (typically the internet), so that remote 
users can operate the software, and store their resulting data either locally on the server or remotely 
on their own or another machine. 

Overall, the software will include one or both of a geomagnetic transmission modeling 
module and a particle flux modeling module. Optionally, the software includes both of these 
modules, and also optionally the software includes one or more of several optional modules, 
including a shielding transport modeling module, a linear energy transfer modeling module, a proton- 
induced single event upset (SEU, a particular type of SEE, but unless noted otherwise used 
interchangeably herein with SEE) rate modeling module, and a direct ionization induced single event 
upset rate modeling module. 

Examples of each of these modules have been produced, and will be discussed in greater 
detail below. The exemplary module for modeling geomagnetic transmission is referred to herein 
as GTRN. The exemplary module for modeling particle flux is referred to herein as FLUX. The 
exemplary module for modeling transport of particles through solid shielding is referred to herein 
as TRANS. The exemplary module for modeling linear energy transfer to a target structure (typically 
an integrated circuit) is referred to herein as LETSPEC. The exemplary module for modeling proton- 
induced SEU event upset rates in specified devices is referred to herein as PUP. The exemplary 
module for modeling SEU rates induced by direct ionization is referred to herein as HUP. 
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Exemplary FORTRAN code for these exemplary modules appears infra. Together, the software for 

these exemplary modules is referred to herein as CREME96. 

Referring to FIG. 1 , one sees that CREME96 produces two types of outputs: proton-induced 

SEU rates, and direct ionization induced SEU rates. Each of these will represent an upset rate for 

a given electronic device in a satellite with a given shielding, under given environmental conditions, 

for a satellite in a given orbit or orbit segment. 

Referring to 100, in the case where the specified satellite orbit is at or above a given altitude, 

nominally geosynchronous, the software will evaluate the fluxes of SEU-inducing particles at the 

outer surface of the satellite, from a number of sources including galactic cosmic rays (GCR), 
anomalous cosmic rays (ACR), and solar heavy ions. The environment at such satellite orbits is 
referred to herein as the "non-trapped energetic particle environment", since it is essentially 
equivalent to the environment at other locations that are at about 1 AU from the sun, but outside of 
the volume where trapped or quasi-trapped particles would be found in significant numbers. 

Referring to 200, in the case where the specified satellite orbit is below a given altitude, 
nominally geosynchronous, before the flux at the outer surface of the satellite is evaluated, the 
shielding effect of the earth's magnetosphere is evaluated. As persons of ordinary skill in the art will 
recognize, the earth's magnetic field will shield out a significant fraction of charged particles. Thus, 
it is necessary to first subtract from the flux in the non-trapped energetic particle environment that 
fraction of the flux that will be shielded by the magnetosphere. 
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Referring to 300, the external hull of a satellite (typically made from thin aluminum panels) 
will shield the electronics of the satellite from a fraction of SEU-inducing particles. Thus, deter- 
mining an SEU rate for microelectronics on a given satellite will require taking this shielding into 
account by correcting the flux for the energy lost in said shielding and for the fragmentation of the 
incident particles within said shielding, and subtracting from the flux of SEU-inducing particles the 
subset of SEU-inducing particles with energies that are not sufficient to penetrate through the hull 
of the satellite. This module of the software performs that function. Although, as indicated by the 
dashed lines in the flowchart, this module may be bypassed, it is highly preferred to not bypass this 
module. 

Referring to 400, although an approximation of the satellite hull shielding may be obtained 
by assuming that the shielding is uniformly distributed about the electronics of interest, more 
realistic results may be obtained by specifying the particular distribution of shielding about the 
electronics for a particular satellite. This module of the software performs that function. 

Referring to 500, after the software measures the transmi ssion of SEU-inducing particles, the 
software will measure the effect of SEU-inducing protons on the electronics of interest. 

Referring to 600, after the software measures the transmission of SEU-inducing particles, the 
software will calculate the linear energy transfer to the electronics of interest. This module of the 
software takes an input file in the form of a particle flux vs. kinetic energy table for each elemental 
species, and generates an output file in the form of a function specifying the combined flux of all 
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these species as a function of the rate of energy deposition. This is an intermediate result for the 
calculation of SEU rates attributable to direct ionization. 

Referring to 700, this program module takes the output linear energy transfer spectrum from 
600, combines this with parameters for the particular electronic device, calculates a SEU rate for the 
device attributable to direct ionization. The parameters are typically based on ground tests of the 
sensitivity to a particular chip to SEU-inducing radiation. 

Having given an overview of this preferred embodiment of the invention, a more detailed 
description of the particular software modules making up this preferred embodiment is given. 

Geomagnetic Transmission Function 

Although it has been known that the magnetosphere shields cosmic rays, quantifying this 
effect to a degree of accuracy that could lead to accurate models of the effects of cosmic rays on 
electronics within the magnetosphere has proven elusive. One difficulty has been the lack of 
accurate models for the magnetosphere. Another difficulty has been that even if an accurate model 
of the magnetosphere were available, modeling the shielding effect of the magnetosphere requires 
an inordinate amount of computer processing time. 

Referring to FIG. 2, this figure shows a representative path of a cosmic ray through the 
earth's magnetosphere (not shown). This path is highly irregular, including numerous changes of 
direction, and covering widely-spaced areas over earth's surface. Determining whether a particle 
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with such a convoluted path will intersect a body as small as a typical satellite is a considerable 
challenge. 

B, GRF is the International Geomagnetic Reference Field promulgated by the International 
Association of Geomagnetism and Aeronomy. It is the internationally accepted standard expression 
for the earth's internal magnetic field. It has been discovered that the B IGRF is insufficient, and that 
the most well-known model of the magnetospheric B fields, the model provided by Tsyganenko et 
al., "A Magnetospheric Magnetic Field Model with a Warped Tail Current Sheet", Planet. Space 
Sci. ,37(1) 5-20 (1 989), is likewise insufficient. Corrections to this model are included in the present 
invention. 

It has been discovered that within the earth's magnetosphere, the total magnetic field can be 
modeled by; 

B = B IGRF (r,f) + B TSYG (Kp,r,f) 
where B IGRF is the International Geomagnetic Reference Field promulgated by the International 
Association of Geomagnetism and Aeronomy, and B TSYG is the modified Tsyganenko field given by 
the sum of B^) + B Ymg < T > + B Zmg ™ + B^ + B Ymg < RC > + B Zmg < RC > + Bxsm < c > + B Ysm < c > + Bzsn / C > + 
B xsm (M) + B Ysm (M) + Bzm m > where coordinates in the solar magnetospheric system are denoted sm, 
and where coordinates in the solar magnetic coordinate system are denoted mg; and where 

Bxsm^ = Qt Xsm 

B Ysm m = Q t y sinZr ; 
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z s (x,y,^)= 0.5 tan ^(x+ R c - -J(x+ R c ) 2 + 16 

-Gsin^-y 4 (y 4 + L 4 y )"', 
where S T RC = ^p 2 + (a TRC + ^ RC f , 

^T.RC = V 2 ^ + ^T,RC ' 

D T = D 0 + Sy 2 + y T hr{x) + 



where C, = -98.72 when Kp = 0, 0^, C, = -35.64 when Kp = r,l,l + , C, = -77.45 
when Kp = 2 ,2,2 + , C,= -70.12 when Kp = 3',3,3\ 0, 0 + , C, = - 162.5 when Kp = 
4\4,4\ C, = - 128.4 when Kp ;> 5", 
5 where C 2 = - 10014 when Kp = 0,0",C 2 = -12800 when Kp=r,l,l + ,C 2 = -14588 

when Kp = 2",2,2 + , C 2 = - 16125 when Kp = 3-,3,3 + , 0, 0 + , C 2 = - 15806 when Kp = 
4\4,4\ C 2 = - 161 84 when Kp * 5", 

where a T = 13.55 when Kp = 0, 0\ a T = 13.81 when Kp = 1\1, 1\ a T = 15.08 when 
Kp = 2-,2,2 + , a T = 15.63 when Kp = 3\3,3 + , 0, 0 + , a T = 16.1 1 when Kp = 4-,4,4 + , a T 
10 = 15.85 when Kp * 5", 
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where D 0 = 2.08 when Kp = 0, 0 + , D 0 = 1.664 when Kp = 1,1,1*, D 0 = 1.541 when 

Kp = 2\2,2 + , D 0 = 0.9351 when Kp = 3",3,3 + , 0, 0 + , D 0 = 0.7677 when Kp = 4',4,4 + , 

D 0 = 0.3325 when Kp * 5", 

where R,. = 8 R E and L y = 10 R E , 

Bx^QrcXZ,; 
B Y (RC) = QrcYZ,; 

n (RO r 2(a RC + 4 C ) 2 - p 2 dz, Rr <?z« <?D RC 

Bz - Cs s RC + B * 77 +B - 77" QrcDrcX "7x" ; 



where Q RC = SC^Rc'SRc'CaR^Rc) 
D RC = Do+YRcMxKYih^x) 
h T.Rc = 0-5[l+x(x 2 +) 

h TRC = 0.5[l + x(x 2 + L 2 TRC )- ,/2 ], 

h, = 0.5(1- (x+ 16)[(x+ 16) 2 + 36r ,/2 } 5 

C 5 (Dst) = - 10220+408.5-Dst 
B XYZ (C) = C 3 (F + x>yiZ + F^ +C 4 (F + x , y>z -F-„ yiZ ), where 
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and where L y = 10.0, D x =13.0, L RC =5.0, L T =6.30, y T =4.0, 6=0.010, Yl =1.0, R T =30.0, ^=4.0, 
L xc =50.0, and D yc =20.0. 

Additional parameters relevant to the present invention given in the Tsyganenko reference 
include (from page 12 of Tsyganenko, supra): 
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It should be noted that the model for B XSYG differs from the model given in Tsyganenko, 
supra. In the first place, the author has identified several errata to this model. In the second place, 
the present inventors have modified the model, making C 5 a linear function of Dst for all geomag- 
netic activity levels. See Boberg et al. 

Referring to FIG. 3, the Lorentz equation using a corrected model for the geomagnetic B field 
is used to calculate the geomagnetic transmission function, using the method depicted in this 
flowchart. Given the large computational requirements for this calculation, it is preferred to make 
this calculation in advance for as many orbits as desired, and storing the results as a series of 
geomagnetic transmission function tables, with the geomagnetic transmission probability given for 
each particle rigidity (rigidity = momentum/charge), where each table corresponds to the geomag- 
netic transmission function for a particular orbit. The user may then select from several predefined 

December 31, 1997 (1:24PM) 17 



Navy Case No. 78,824 PATENT APPLICATION 

Inventor(s): Adams et al. 

orbits having predefined geomagnetic transmission functions. Alternatively, especially as computer 
processing power increases, in another embodiment of the invention, this calculation of the 
geomagnetic transmission function may be performed after a user has specified an orbit. 

Referring to 300, the calculation of the geomagnetic transmission function will typically 
begin by specifying a particle rigidity in GV, a range of particle arrival directions by specifying a 
range of azimuth angles between 8 max and 6 min and a range of zenith angles between d> and <b • 

• max t mm' 

specifying an earth orbit, and specifying a geomagnetic activity level by specifying values for Kp and 
Dst, where Kp = [IATME or IAGA] planetary Kp index and Dst = the hourly equatorial Dst index. 
A standard reference for the definition of Kp is Bartels, J., "The standard index, Ks, and the 
planetary index, Kp", IATME Bulletin 12b, 97 IUGG Pub. Office, Paris, 1949. IATME is the 
International Association of Terrestrial Magnetism and Electricity. IAGA is the International 
Association of Geomagnetism and Aeronomy. IUGG is the International Union of Geodesy and 
Geophysics. A standard reference for the definition of Dst is Sugiura, M., "Hourly values of 
equatorial Dst for the IGY", Ann. Int. Geophys. Year, 35, 49, 1964. The Dst index was adopted by 
IAGA at the 1969 meeting in Madrid. IAGA Bulletin 27, 1969, 123, resolution 2. All of the 
references cited in this paragraph are incorporated herein by reference, in their entireties, for all 
purposes. 

For each value of particle rigidity, the transmission probability will be calculated, and each 
of these results will be compiled into a table giving the geomagnetic transmission function for a 
given orbit. A separate geomagnetic transmission function table will be prepared for each specified 
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geomagnetic activity level. The user will also specify an arbitrary number of arrival directions per 
orbital step M^. 

Referring to 305, the user will specify an arbitrary number of orbital steps N spo , where N sp0 
is the number of steps in one orbital revolution of the satellite around the Earth's spin axis. The user 
will also specify a required accuracy for the numerical integration. This embodiment of the 
invention relies on tracing the particle trajectories backwards from a discrete number of evenly 
spaced satellite orbital positions. Persons of ordinary skill in the art will recognize, however, that 
alternative embodiments might be employed, such as embodiments using a Monte Carlo simulation, 
to trace particle trajectories backwards from randomly selected satellite orbital positions. 

Referring to 310, the user will specify the satellite orbit. 

Referring to 315, the software will calculate the orbital period from the orbit parameters, 
using standard orbit generation routines. See Adams et al., NRL Memorandum Report 5099, supra. 
The software will also calculate a maximum number of steps, and a time per step for the satellite to 
travel through a step. 

Referring to 320, the software will calculate the position of the satellite after each step. See 
Adams et al., NRL Memorandum Report 5099, supra. The software will also begin tracing the 
position of a particle incident on the satellite by initializing the particle trajectory time. 

Referring to 325, for each of the specified number of arrival directions per orbital step, the 
software will select a random particle arrival direction. Persons of ordinary skill in the art will 
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recognize, however, that alternative embodiments might be employed, such as embodiments where 
a number of regularly spaced arrival directions are used. 

Referring to 330, the software will calculate the particle's final position and velocity vector 
v from the satellite's latitude, longitude, altitude, and the particle's rigidity, final 6, and final <t>. 

Referring to 330, 340, 345, the software will calculate a new particle position by performing 
a numerical integration of the Lorentz equation. Although this preferred embodiment uses the 
Bulirsch-Stoer numerical integration technique, persons of ordinary skill in the art will recognize that 
other numerical integration techniques are available. This numerical integration is performed on a 
modified Lorentz equation, where Q is replaced by -Q and f is replaced by -f. As persons of 
ordinary skill in the art will note, the Lorentz equation is invariant if one makes this substitution. 
This significantly simplifies the numerical integration, because the solution in this case will represent 
tracing a particle backwards in time from the satellite to the boundary of the magnetosphere (for 
allowed trajectories) or earth's surface or atmosphere (for forbidden trajectories). This will typically 
entail far fewer calculations than if particles incident on the magnetosphere were traced forward in 
time to see if they intersected the satellite. However, persons of ordinary skill in the art will 
recognize that this tracing forward in time is an equally valid method. 

Referring to 350, 355, 360, the software will determine whether, after this iteration of the 
numerical integration, the particle's new position has placed it within the atmosphere or the solid 
earth (representing a forbidden trajectory), beyond the magnetosphere as given by the model 
(representing an allowed trajectory), or still within the magnetosphere (representing an undetermined 
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trajectory). If the particle is still within the magnetosphere, another iteration of the numerical 
integration is performed, and a new particle position determined. If the particle is not still in the 
magnetosphere, i.e., particle has reached either earth's surface or atmosphere (forbidden) or the 
outside of the magnetosphere (allowed), the software will store this result. This preferred 
embodiment of the software has an error trapping technique, where if the number of particle steps 
becomes very large, the particle is forced to an allowed or forbidden transition. 

Referring to 365, if the software has reached a solution for a given particle arrival direction 
at a given orbital step, the software will move on to trace the particle arriving at the next arrival 
direction for that orbital step, if any. 

Referring to 370, if the software has reached solutions for all the arrival directions for the 
satellite at a given orbital step, the software moves the satellite along to the next orbital step, if any. 

Referring to 375, after all particle tracings for all orbital steps have been performed, the 
results are compiled as a geomagnetic transmission function of the form GT(R) = NJN max , where 
GT is the probability of a particle being transmitted through the magnetosphere, R is particle rigidity, 
and N tt and N max are a number of particles transmitted and a total number of particles, respectively. 
Typically, this function will be embodied in a lookup table. 

Calculating a geomagnetic transmission function for a given orbit according to the above 
method is an intensive process. Computing geomagnetic transmission functions for two particular 
orbits (the proposed space station orbit and a typical space shuttle orbit) required several days of 
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processor time on a RISC workstation. Shorter computation cycles could be made on a super- 
computer, but supercomputer time is in chronic short supply. Accordingly, it is preferred to calculate 
geomagnetic transmission functions for preselected orbits in advance, and to store these geomagnetic 
transmission functions for later use by users. This is the approach taken by the inventors in their 
CREME96 software. In cases where a user specifies an orbit that is not among the precalculated 
orbits, the CREME96 software will use an earlier model that uses geomagnetic cutoff tables (see 
NRL Memorandum Report 5099, supra) rather that geomagnetic transmission functions, and warns 
the user of this substitution. 

Alternatively, a user would specify an orbit, and if the software did not have a precalculated 
geomagnetic transmission function for this orbit, the software would calculate this geomagnetic 
transmission function. This embodiment is anticipated to become increasingly desirable as 
processing capacity increases. 

Referring to FIGS 4, 5, 6, and 7, these plots show the geomagnetic transmission function for 
a common space shuttle orbit (circular orbit, 450 km altitude, 28.5° inclination) under "quiet" 
(FIG. 4) and "stormy" (FIG. 5) conditions, and the geomagnetic transmission function for the 
proposed space station orbit (circular orbit, 450 km altitude, 51.6° inclination) under "quiet" (FIG. 6) 
and "stormy" (FIG. 7) conditions. "Quiet" magrietospheric conditions were taken herein to be 
Kp = 2 and Dst = -15 nT, while "stormy" magnetospheric conditions were taken herein to be Kp * 5 
and Dst = -300nT. One sees that (1 ) transmission drops off sharply with decreasing rigidity, and (2) 
the differences between stormy and quiet conditions are seen most at lower rigidities. 
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The embodiment of the invention implemented in the CREME 96 software operates by 
prompting a user to specify an orbit and magnetospheric conditions, and then looking up the 
appropriate geomagnetic transmission function table for the orbit and magnetospheric conditions. 
The CREME96 software includes lookup tables corresponding to each of FIGS. 4 through 7. This 
table is then applied to a flux of SEU-inducing particles in the non-trapped energetic particle 
environment (calculated by another part of the CREME 96 software), to determine the flux of such 
SEU-inducing particles at the outer surface of the satellite. In other words, the flux from the non- 
trapped energetic particle environment will be multiplied by the probability (as embodied in the 
geomagnetic transmission function) that the particles will be transmitted through the magnetosphere, 
to get the flux at the outer surface of the satellite. 



Particle Flux Modeling 

Four principal classes of particles are considered to be responsible for most SEUs: galactic 
cosmic rays (GCRs), anomalous cosmic rays (ACRs), solar protons, and solar heavy ions. See NRL 
Memorandum Report 4506, supra. 

To estimate the effect of these classes of particles on orbiting microelectronics, it is first 
necessary to estimate the exposure of the microelectronics to these particles. Thus, a feature of any 
complete model for the effect of cosmic rays on microelectronics will be the modeling of the flux 
of each of these types of particles. 

In the present invention, particle fluxes are determined at the near earth environment. As 
used herein, the near earth environment refers to the environment found at locations approximately 
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1 AU from the sun, but outside of earth's magnetosphere or other localized phenomena that would 
perturb particle fluxes. After the particle fluxes for each type of particle is determined by the FLUX 
module of a preferred embodiment of the present invention, other modules may be used to model 
how the charged particles will induce single event upsets, including what shielding effects (including 
shielding by the magnetosphere and shielding by the hull of the spacecraft) will reduce this flux at 
the component level. 

In a preferred embodiment of the invention, the flux of GCRs is modeled by the galactic 
cosmic ray model described by Nymmik et al. See "A Model of Galactic Cosmic Ray Fluxes", Nucl. 
Tracks Radiat. Meas. 20 427, Nymmik et al. (1992), incorporated herein by reference, in its entirety, 
for all purposes. Although not essential to the present invention, the authors' updated version of this 
model is described at http://www.npi.msu.su/gcrf/standart/ISO_WD_15390.html, incorporated 
herein by reference, in its entirety, for all purposes. 

In a preferred embodiment of the invention, the flux of ACRs and solar protons are modeled 
by the model described in Tylka et al, supra. 

In a preferred embodiment of the invention, the flux of solar heavy ions is modeled as 
depicted in FIGS. 8 and 9. This routine depicts how a solar heavy ion flux for a given atomic species 
at a given kinetic energy, under conditions that are selected from one of three available models: 
worst day, worst week, and peak flux. As used herein, the worst week solar heavy ion model refers 
to the model based on particle fluences averaged over a 1 80 hour span of the solar event beginning 
at 1300 UT on October 19, 1989. As used herein the worst day solar heavy ion model refers to the 
model based on particle fluences averaged over an 1 8 hour span (beginning at 1300 UT on October 
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20, 1989). As used herein, the peak flux solar heavy ion model refers to the model based on the 
highest fluxes averaged over 5 minute intervals reported by the Geostationary-orbiting Operational 
Environmental Satellite (GOES) satellite. 

Referring to 800, this routine takes as inputs the atomic number of the selected species, and 
the kinetic energy for that species whose flux is to be determined. The routine also takes as inputs 
the user selected baseline model (worst week, worst day, or peak flux). 

Referring to 805, if the selected species has atomic number greater than 20, iron is selected 
as the model spectrum for the flux of this species 810. Otherwise, 815, oxygen is selected as the 
model spectrum for the flux of this species. 

Referring to 820, the next step is to look up an elemental breakpoint EB2. This breakpoint 
represents a value for the kinetic energy where the relationship between the flux and the kinetic 
energy changes to a degree where a different mathematical expression for this relationship is 
warranted. Values for EB2 are given in Table S 1 : 



Table SI : Values for Elemental Breakpoint EB2, in MeV/nuc 



^^■■■^^^ Baseline model 
Elemental 

spectrum model ' — 


worst day or peak flux 


worst week 


iron 


24.23 


19.90 


oxygen 


15.94 


12.89 



Optionally, this lookup may be performed later in the routine, after the IF step at 825. 



Referring to 825, a special case exists for when the baseline model is the worst week model, 
the model element is Fe, and the kinetic energy being evaluated is above a certain cutoff point 
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E sp = 1 27.93 MeV/nuc. In this case, referring now to FIG. 9, values for A, p and G 5p are looked up from 
Table S2 and Table S3, respectively, and the unsealed flux is calculated to be A sp xEn" 0sp , 900,910. 



Table S2: Values for A. n : 



* — Baseline model 
Elemental ^^■'^^^ 
spectrum model — 


aw 

worst day or peak flux 


worst week 


iron 


not applicable 


3.16814xl0 6 


oxygen 


not applicable 


not applicable 


Table S3: Values for G, n : 


' ^^^^ Baseline model 
Elemental ^" — 
spectrum model — 


worst day or peak flux 


worst week 


iron 


not applicable 


2.861 


oxygen 


not applicable 


not applicable 


Referring back to FIG. 8, if En > EB2 (835), values for A 3 and y si are retrieved from tables 
S4 and S5 ? respectively (840), and the unsealed flux is calculated to be A 3 xEn" YSi (845). 

Table S4: Values for A, 


^^^-^^ Baseline model 
Elemental * N,,, * > --^ 1 ^^ 
spectrum model — 


worst day or peak flux 


worst week 


iron 


0.252948* 10'° 


0.249719* 10 9 


oxygen 


0. 106702* 10 10 


0.667628* 10 9 
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15 



i« 20 

: 1 



^^^^^ Baseline model 
Elemental ^^^^ 
spectrum model ^^^^ 


* *' 

worst day or peak flux 


worst week 


iron 


4.52970 


3.7610 


oxygen 


4.14060 


3.76850 



Otherwise, for En*EB2, values for A 2 and G sj are retrieved from tables S6 and S7, 
respectively (855), and the unsealed flux is calculated to be A 2 exp(-GxEn l/4 )xEn 1M (860). 



Baseline model 
Elemental ^-^^^ 
spectrum model ***^ <iii ^^ 


i ■ — L 

worst day or peak flux 


worst week 


iron 


1.8991xl0 8 


3.0372xl0 8 


oxygen 


4.9518wxl0 8 


1.1307xl0 9 


Table S7: Values for G c; 


Baseline model 

Elemental 

spectrum model ---^^ 


worst day or peak flux 


worst week 


iron 


5.70 


5.70 


oxygen 


5.70 


5.70 



For heavier species (IZ*3), the flux must be scaled by looking up scale factors for the species being 
25 analyzed and the model element (if different from the species being analyzed). These scale factors 
are taken from Table S8: 
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Table S8: Scale factors SF for atomic species with IZ^3 





Atomic 
Number 


Scale Factor 


Atomic 
Number 


Scale Fac- 
tor 


Atomic 
Number 


Scale Fac- 
tor 


Atomic 
Number 


Scale Fac- 
tor 


Atomic 
Number 


Scale Fac- 
tor 


Atomic 
Number 


Scale Fac- 




5 


0 


22 


4.377 
xlO' J 


39 


4.878 
xlO** 


56 


4.878 

xlO* 6 


73 


2.195 
xlO" 8 


90 


4.878 
xlCr* 


5 


6 


4.704 

xlO"' 


23 


4.088 


40 


1.22 
xlO* 5 


57 


4.878 
xlO* 7 


74 


2.439 
xlO' 7 


91 


0 




7 


1.2059 
xio' 


24 


1.65 
xlO' 1 


41 


9.756 
xlO* 7 


58 


1.22 
xlO* 6 


75 


4.878 
xlO"* 


92 


2.927 
xio 4 




8 


1 


25 


5.625 
xio 3 


42 


4.878 
xio* 6 


59 


1.951 
xio* 7 


76 


7.317 
xlO* 7 


3 


0 




9 


4.560976 
xio- 3 


26 


1 


43 


0 


60 


9.756 
xlO' 7 


77 


7.317 
xlO* 7 


4 


0 




10 


2.1312 

xJO*' 


27 


1.303 
xtO' 2 


44 


2.195 
xlO" 6 


61 


0 


78 


1.463 
xio* 6 






10 


11 


1.744715 . 

xlO" 2 


28 


3.172 
xlO* 2 


45 


4.878 
xio 7 


62 


2.439 
xlO* 7 


79 


2.439 
x 10' 7 








12 


2.0624 


29 


3.048 
xlO" 4 


46 


1.463 
xlO" 6 


63 


9.756 

xlO"* 


80 


2.439 
x in* 7 








13 


1.826829 


30 


7.457 
xlO* 4 


47 


4.878 
xio 7 


64 


4.878 

xlO' 7 


81 


2.195 
x in- 7 








14 


3.5935 

xlO' 1 


31 


4.878 
xlO* 3 


48 


1.707 
xio- 6 


65 


7.317 

xlO"* 


82 


2.439 
x 10" 6 








15 


2.279675 
xlO" 4 


32 


1.22 
xlO" 4 


49 


2.195 

xlO 7 


66 


4.878 
xlO* 7 


83 


1.463 
x 10' 7 






15 


16 


9.758 

xlO" 2 


33 


7.317 
xio- 6 


50 


4.878 
xlO* 6 


67 


9.756 
xlO 4 


84 


0 








17 


1.680488 
xlO- 4 


34 


7.317 
xlO* 5 


51 


3.415 
xlO* 7 


68 


2.439 
xlO' 7 


85 


0 








18 


1.771545 
xlO* 3 


35 


9.756 
xlO* 6 


52 


7.317 
x 10* 6 


69 


4.878 
xlO** 


86 


0 








19 


3.644715 
xlO-* 


36 


4.878 
xlO 3 


53 


1.463 
xio- 6 


70 


1.951 
xio 7 


87 


0 








20 


4.826 
xio* 2 


37 


7.317 

xlO" 6 


54 


6.585 
xlO* 


71 


4.878 
xlO 4 


88 


0 






20 


21 


2.929 
xlO" 4 


38 


2.439 
xio 5 


55 


4.878 
xlO" 7 


72 


1.951 
xlO" 7 


89 


0 







As embodied in preferred CREME96 software according to the invention, this routine is run 
for each atomic species in the solar heavy ion flux, and for each value of kinetic energy in the solar 
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heavy ion flux. To determine the total solar heavy ion flux, the fluxes for each kinetic energy of each 
species are added together. 

Internet Implementation 

The preferred CREME96 software embodiment of the present invention is large and 
complex. Moreover, it is anticipated that as models used in the software are improved and/or 
updated, it will be necessary to perform maintenance on this software. Accordingly, as noted supra, 
a preferred embodiment of the invention includes a network interface for allowing remote users to 
access and run the software. Since the internet is a stateless system, a user does not maintain a 
connection to the server as users in a client-server environment would. One of the goals of the 
internet implementation of the present invention was to provide users with as much of the 
functionality of a client-server based system as possible in the stateless system of the internet. 

The preferred system for allowing a large number of users to remotely access and run 
software such as CREME96 is a system connected to the internet, or equivalently a system connected 
to an intranet. Accordingly, the preferred embodiment of the invention is a server computer (or 
equivalently, a plurality of server computers connected in parallel) with the CREME96 software, 
with PERL script software (or equivalently Java scripts, or scripts written in some other script 
language) for controlling input to and output from the CREME96 software, as well as other 
functions, and with listener software and other hardware and software necessary to permit remote 
users with internet access to access the CREME96 software through their respective web browsers. 
Listings for the preferred software according to the invention are found infra, 
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Referring to FIG. 10, this figure shows a main menu for the preferred CREME96 software 
according to the invention. This menu may be accessed through the CREME96 web site maintained 
at the Naval Research Laboratory having a URL http://crsp3.nrl.navy.mil/creme96/. This page is 
generated by a PERL script in response to a registered user transmitting to the server a valid 
username and password. The PERL script will compare the transmitted username and password to 
a list stored on the server, and if a matching entry is found, the script will generate a CREME96 main 
menu page customized for that user. 

One of the advantages of CREME96 over other web-based applications is that it stores 
information for each user on the server. This obviates the need for cookies or other downloads to 
the user's computer. Many users find the downloading of cookies to be objectionable. Moreover, 
keeping this information on the server may speed up overall operation, because there is no need for 
the server to download this information to the user, and the user to subsequently upload it to the 
server. Thus, when a user enters a valid username and password, the script will search a section of 
the hard drive set aside for that user for any User Request Files (URFs). The user will then be able 
to pick any of the user's previously stored URFs to use as a data file for a selected routine. 

The user has the option, from this main menu, of either running one or more modeling 
routines, using a selected URF, editing or creating an URF, or running utilities. 

The available routines include the GTRN routine for calculating a geomagnetic transmission 
function as described supra, the FLUX routine for calculating a flux of SEU-inducing particles, as 
described supra, the TRANS routine for calculating the transmission of SEU-inducing particles 
through solid shielding, as described supra, the LETSPEC routine for calculating the linear energy 

December 31, 1997 (1:24PM) 30 



Navy Case No. 78,824 PATENT APPLICATION 

Inventor(s): Adams et al. 

transfer of charged particles to electronic devices, as described supra, and the PUP and HUP routines 
for calculating the SEU rates for electronic devices exposed to SEU-inducing particles. 

The user will select routines to run by clicking a check box next to each of the routines to be 
run, and picking an URF for each routine to be run from the list of URFs that user has previously 
stored for that routine. When the user sends this message to the server, by pressing the GO button, 
a script will generate a system command to run each of these selected routines, using the selected 
URF as an input file for that routine. As persons of ordinary skill in the art will recognize, a system 
command is a computer command at the system level. In this case, the system command is to run 
a program with specified inputs. 

Alternatively, if the user selects a radio button next to one of these routines ancfclicks the GO 
button, this message from the user will trigger a script to generate and transmit to the user an HTML 
page with inputs for the data fields in the corresponding URF. If the user has first specified an 
existing (i.e., stored on the server) URF, the values for the fields in this URF will be included in the 
generated HTML page, allowing the user to edit the values for the fields in the URF. The URF file 
name is included in the URF fields, thereby providing the convenience of creating a new URF 
derived from editing a previous URF, and allowing file storage of both the original and modified 
URF. 

One of the features of the preferred embodiment of the invention is the connectivity between 
the routines. Referring back to FIG. 1, one sees that the modules of the software are typically 
executed in the indicated order. In the preferred embodiment, the modules receive input in the form 
of a URF. The modules also receive input from a file generated by the preceding module (this file 
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specified by the URF), and/or send output to a file for use by the subsequent module (this file 
likewise specified by the URF). For instance, the GTRN module will generate a .GT* output file. 
The FLUX module may use this .GT* file as an input, and will output a .FLX output file. The 
TRANS module will typically use a .FLX file as an input, and output a .TFX file. The LETSPEC 
module will typically use a .TFX file (or optionally a .FLX file) as an input, and will output a .LET 
file. The PUP module will typically use a .TFX file (or optionally a .FLX file) as an input, and will 
output a .PUP proton-induced upset report file. The HUP module will use a .LET file as an input, 
and will output a .HUP direct ionization induced upset report file. 

Accordingly, by specifying a URF for each of these routines, selecting each of these routines 
to run, and clicking the GO button, a user will send a message to the server that will cause the scripts 
to run each of these routines in order, using the outputs from the preceding routines and generating 
inputs for the subsequent routines. The result of these calculations will be estimates of the single 
event upset rates for the specified electronics in the specified orbit under the specified conditions. 
The scripts will generate and transmit an HTML page with these results to the user, who may then 
use the utilities to download the results, create plots, or perform other functions. 

Having described the invention, the following examples are given to illustrate specific 

applications of the invention, including the best mode now known to perform the invention. These 

specific examples are not intended to limit the scope of the invention described in this application. 

Example 1: 
CREME96 
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Code for CREME96, a preferred embodiment of the present invention, follows. The software 

may be accessed at its internet site with URL http://crsp3.nrl.navy.mil/creme96/. 

Obviously, many modifications and variations of the present invention are possible in light 
of the above teachings. It is therefore to be understood that, within the scope of the appended claims, 
the invention may be practiced otherwise than as specifically described. 
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